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CLASFY6 PROGRAM 


1 .0 SCOPE 

This document contains the description of the CLASFYG computer program. 
The program calculates coefficients to fit the model 

G(t) ■ A /t_\ “ -e(tj - t^) for t > t„ 

\tj t 0 

and 

6(t) ■ A for t < t^ 

where G Is the greenness and A Is the soil line, to the time trajectory 
of each pixel In a specified region of a segment. A disk file Is written 
with a Universally formatted header record and line records recording the 
Oil 3t and CHI-Square value for each pixel on the scan line. 


1-1 


2.0 APPLICABLE OOCUFCNTS 


The following documents form a part of this specification: AD 63-2457-3308-01 
Transferring of Badhwar Software. 


AD NAS 9-15200 Technical Memorandum Format Specifications for LACIE 
(Phase III) and Accuracy Assessment Computer Data Products. 


Ill 
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3.0 SYSTEM DESCRIPTION 


3.1 HARDWARE DESCRIPTION 


Th« softMare for CLASTYG Is operational on the IBM 3031 computer at Purdue. 


3.2 CLASFY6 PROCESSOR SYSTEM FLOWCHART 


The system level data flow diagram is shown in figure 3.1. 


The following system flowchart assumes the use of the EXEC file described 
in Appendix C. 


.\cqui»ition duta files (S*8) 
on nultifile upe(s) 
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Th« Inputs and outputs shown on Figurt 3.1 are described briefly here. 
A more detailed description Is given In section 3.4. 


Acquisition data files 


User Information file 


Data from RT4E data base 


Sumnary Report 


Coefficients File 


- Header record and 4 channels of data for 
each of 5 to 8 acquisitions (See Section 
3.4.3). 


- User Information file supplying crop, 
Initial guesses for greenness values, 
training field vertices and optional 
test field vertices (See Section 3.4.1). 

file ' Segment number, number of acquisitions, 

Landsat numbers and sun angles (See Section 
3.4.2). 

- Report of date and time of processing, 
crop vertices for training field, 
acquisition dates, all statistical data 
created by CLASFYG and time taken to 
execute the program. 

- File of one header and 117 data records, 
each data record containing the Alpha, 

Bct^, t^ and CHI square values for each- of 
the 196 pixels In a line of acquisition 
data (See Section 3.4.4). 
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CLASFY6 Hierarchy Diagram 


CLASFYG (CLASFYG)* — CROPG 


(CROPG),- 


— MOVSYM 
->NUHBFL 
— NUWR 


(CROPG) 

(CROPG) 

(CROPG) 


— ROEXEC 
— ROFLES 
— RECT 
— SDEV 
— TRJXM 


(ROEXEC: 

(RDFLES)-RDVERT 

(RECT) 

(SDEV) 

(TRJXM)p-F2ERVG 
— F2NCTG 
— FCHISQ 
— KTINVl 


(RDFLES) 


( F2NCTG) 
(F2NCTG) 
( F2NCTG) 
(MTINVl) 


* Parentheses enclose the name of the file which contains the source 
statements for the named program which precedes them. For example, 

PROGRAM (Source statement file name) 


3.3 SOFTWARE DESCRIPTION 


Tht CLASFY6 program fits a modtl to the time trejcctory of each pixel In a 
specified region of a segment (usually this region Is the whole segment). 

The time trajectory consists of a plot of the greenness values for a 
particular pixel on various acquisitions as a function of time. The model 
Is given by 

6(t) - A ® ^ -B(t^^ - t^ ) for t > t^ 

and 

G(t) • A for t <. t^ 

where G Is the greenness and A Is the "soil line". The CLASFYG program 

produces estimates of the parameters A, a. S and t^ for each pixel which 

give the best fit to the trajectory defined by the data for that pixel. 

2 

It also calculates a value for x which describes the goodness of fit. 

The first stage of the processing consists of determining starting values 
for A, a, 6 and t^ which are subsequently used as starting points for 
fitting the model to pixels In the test area. 

These starting values are determined using the data In a training field 
selected by the user and a set of Initial "guesses" for these parameters 
also Input by the user. The training field should be a field of the crop 
of major Interest. Also during this first stage, an array SGMX Is calculated 
which contains the standard deviations of the pixels In the training field 
for each acquisition used. These are used In weighting the acquisitions 
when fitting the model to the pixels In the test area. 

The second stage of the processing consists of fitting the model to each 

pixel in the test area starting the iteration procedure for each pixel 

with the starting values calculated In phase I. When this Is completed 

2 

the final values for a, B. t and x For each pixel are output as a 

0 

formatted disc' file. Also a sunvnary report Is printed on the line printer. 
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3.4 FILE DESCRIPTION 


3.4.1 USER INFORMATION ^ILE (Segment number) USERG A 

The following records describe the name of the crop, the Initial guesses 
for the coefficients for greenness and the training and test fields and 
allows the user to Input any number of free- form comments. 

In record types 1-5 the **keyword" Is In columns 1-4 and any parameter or 
conrent data Is In columns 11 through 72. Numbers In a series are 
separated by comnas; blanks are optional. 

Record types 6 and 7 which describe the fields are blank In columns 1 
through 10 and the field definition Is In columns 11-72. 

Record type 6 must follow the "*EN0*' record and must proceed record type 7. 

Record type 7 Is optional, but If It exists It must follow record type 6. 

1 . "COMMENT" record 

EXAMPLE: COMMENT THIS IS AN EXAMPLE EXECUTION OF THE PROGRAM 

These records allow the user to enter any number of free-form comments. 

2. "CROP" record 
EXAMPLE: CROP CORN 

This record specifies the 4 character name of the crop of Interest. 

The name begins at the first non-blank character. 

3. "GREENNESS" record 

EXAMPLE: GREENNESS 12.00, 11.50, 1.50, 1.50 

This record specifies the Initial guesses for the A, Alpha, Beta and t^ 
values of the crop In the training field. 
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4. "SUNANGLE" rtcord 

EXAMPLE: SUNANGLE 45. 47. 51, 46. 42 

1h1s rtcord sp«c1f1ts the sun angles and Is used only If the sun angles 
In the data base are wrong. 

5. "*EN0" record 

This record specifies the end of the user Information except for the 
definition of the training field and test field. 


6. "TRAINING FIELD" record 

EXAMPLE: (1.1), (97.100), (97.120). (125,110), (100,100) 

This record specifies the coordinates of the corners of the training 
field. The format of the line Is as follows: 


(1. 

1) 

(x^ , 

’l) 



(Xj. 

''3> 

(x^ » 

’'xl 


Dumny variable to be used for sample skip factor and 
line skip factor; not currently used. 

Unper leftmost corner of training field expressed as 
(sample, line). 

Upper rightmost corner of training field expressed as 
(sample, line). 

Lower rightmost corner of training field expressed as 
(sample, line). 

Lower leftmost corner of training field expressed as 
(sample. 11n«). 


7. "TEST FIELD" record 

EXAMPLE (1,1), (1.1), (196.1). (196,117), (1,117) 

This optional record specifies the coordinates of the corners of the 
test field. The default Is the whole segment. The format Is the same 
as the format for the training field record. 
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.2 AcquisiTiori uata files 

! fulloMing records of the 5 to 8 acquisitions data files are Input to the 
\SFY6 program and are read by RDFLES. 


ese files each contain 4 channels of LANOSAT data for one acquisition, 
an acquisition Is used twice the corresponding file will be duplicated, 
te that these files are all In universal format. Record 1 Is read with 
17{180A1) format statement; records 2-118 are read with a 5(180A1) format 
tatement. 


ecord Bytes 

1-3060 

’ 1-72 

73-268 

269-464 

465-660 

661-856 

3 1-856 


Contents 

Header record In universal format. 

Filler. 

Gray level values for the 196 samples of and 1 
for 1 ine 1 . 

Gray level values for the 196 samples of band 2 
for line 1 . 

Gray level values for the 196 samples of band 3. 
for line 1 . 

Gray level values for the 196 samples of band 4 
for 1 Ine 1 . 

Same data format as record 2; gray level values 
for bands 1-4 for line 2. 


N 


1-856 


Same data format as record 2; gray level values 
for bands 1-4 for line N-1. 


118 


1-856 


Same data format as record 2; gray level values 
for bands 1-4 for line 117. 
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3.4.3 MTA FROH RT&E DATA BASE 


LANOSUN FILE 0 


The following records describe the segment number, the number of acquisitions, 
and the date, sun angle and Landsat numbers for each acquisition. 


Record 

Bytes 

Format 

Contents 

1 

1-4 

A4 

Segment number. 

2 

1-4 

14 

Number of acquisitions 

3 

1-56 

8(2X,A5) 

Acquisition dates. 

4 

1-32 

814 

Sun angles. 

5 

1-32 

814 . 

Landsat numbers. 
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3.4.4 COEFFICIENTS FILE (Segment number) COEF A 

The following records are output from the CLASFYG program: 
Format 


Record 

Bytes 

Format 

Contents 

1 

1-3060 

17(180A1) 

Header record In Universal format. 

2-118 

3136 

18(45A4) 

Coefficients Alpha, Beta, t and CHI-square 


for the 196 pixels on 1 line. 


Detailed format 
Record Words 

1 1-765 

2 1 

2 2 

3 

4 


Contents 

Header record In Universal format. 

Alpha value for pixel 1 In line 1 (computational Integer). 

Beta value for pixel 1 In line 1 (computational Integer). 

t^ value for pixel 1 In line i (computational Integer). 

CHI-square value for pixel 1 In line 1 (computational 
Integer) . 


1+N*4 Alpha value for pixel N+1 In line 1 (computational Integer). 

2+N*4 Beta value for pixel N+1 In line 1 (computational Integer). 

3+N*4 t^ value for pixel N+1 In line 1 (computational Integer). 

4+N*4 CHI-square value for pixel N+1 In line 1 (computational 

Integer) . 


2 


781 Alpha value for pixel 196 In line 1 (computational Integer). 

782 Beta value for pixel 196 In line 1 (computational Integer). 

783 Lq value for pixel 196 In line 1 (computational Integer). 

784 CHI-square value for pixel 196 In line 1 (computational 
Integer) . 
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Cotfficitnts file (continued) 


R ecord 

J 

m 

118 


Words Contents 

1-784 Same data format as record 2; Alpha, Beta, t and 

CHK-square values for 196 pixels of line 2. ° 


1-784 Same data format as record 2; Alpha, Beta, t and 

CHI-square values for 196 pixels of line M. ° 


1-784 Same data format as record 2; Alpha, Beta, t and 

CHI-square values for 196 pixels of line 117? 
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3.5 DETAILED SOFTWARE DESCRIPTION 
3.5.1 CLASFYG 

Purpose 

CLASFYG uses a user determined training field In a full scene to determine 
trajectory coefficients for each band of radiance values to fit the model 
defined In section 1.0. These coefficients Alpha, Beta, t^ and CHI-square 
are written for each pixel to the output file described In section 3.4.3. 

Linkages 

CLASFYG calls CROPG, RDEXEC, CPTIfC, RECT, SDEV, TRJXM and ROFLES. 

Interface 
Calling sequence: 

Not applicable (an EXEC that can be used to call CLASFYG Is described 
In Section 4.0} . 

Calling sequence parameters: 

Not applicable. 

Labeled COMMON parameters: 

Full description of labeled COMMON blocks are contained In Appendix A. 


Element Input/ 


Label 

Variable 

Posi tion 

Output 

Description 

/LINES/ 

LINENO 

2 

0 

Line number. 

/LINES/ 

IDATA(8S6,8) 

3 

I 

Radiance values for 
a maximum of 8 
channels . 

/LINES/ 

NOACOS 

4 

0 

Number of acquisitions 

/MI SC/ 

NOACQ 

1 

I 

Number of acquisitions 

/MI SC/ 

ACQ0T(8) 

2 

I 

Acquisition dates. 

/MI SC/ 

KSEGM 

3 

I 

Segment number. 

/USER/ 

ICR0P(4) 

1 

I 

Name of crop. 

/USER/ 

ZX(4) 

2 

I 

Initial guesses for 
coefficients A, Alpha, 
Beta and t^. 

/USER/ 

SANGGG(8) 

3 

0 

Sun angle correction 
factor. 
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Blank COMMON parameters: 
None. 

Inputs 


Unit 

Type 

Description 

7 

Seq. data 

Data from RT&E data base (See Section 3.4.3). 

11 

Seq. data 

Acquisition data (Sec Section 3.4.2). 

12 

Seq. data 

Acquisition data (See Section 3.4.2). 

13 

Seq. data 

Acquisition data (See Section 3.4.2). 

14 

Seq. data 

Acquisition data (See Section 3.4.2). 

15 

Seq. data 

Acquisition data (See Section 3.4.2). 

16 

Seq. data 

Acquisition datd (See Section 3.4.2) optional. 

17 

Seq. data 

Acquisition data (See Section 3.4.2) optional. 

18 

Seq. data 

Acquisition data (See Section 3.4.2) optional. 

21 

Seq.' data 

User defined parameters (See Section 3.4.1). 

30 

Seq. data 

Scratch file to reformat data. 

Outputs 



Unit 

Type 

Description 

6 

Report 

Summary report . 

8 

Seq. data 

Header and coefficients Alpha, Beta, t and 
CHI-square. ° 

30 

Seq. data 

Scratch file to reformat data. 


Storage requirement 
Not applicable. 


Description 

CLASFYG calls CROPG to read the user specified parameters crop of Interest 
and Initial guesses for greenness coefficients A, Alpha, Beta, t^. 

CLASPTG calls RDFLES to read the header record and the parameters for the 
training field. 

CLASFYG then calls CPTIME to get the current time and date and writes the 
report heading. 
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CLASFYG calculates the mlnlnum and maximum boundaries for pixels In the test 
field, calculates days elapsed from the beginning of the year as a multiple 
of 100 for each acquisition date and sets a switch If an acquisition date Is 
used twice. 

The rectangular field bounded by the minimum and maximum pixel numbers and 
line numbers Is then processed using the following steps: 

A line of radiance values Is read. 

Each pixel In the field Is checked to be In the training field. 

The values for the pixels In the training fields are used to calculate a 
greenness value which Is saved In a matrix. 

When sufficient pixels have been accumulated, the mean and standard deviation 
are calculated for the test field. 

CLASFYG calls TRJXN to calculate Improved values for A, Alpha, Beta and t^ and 
a CHI-square value for the training field. 

The header record Is constructed from constants In the program. Input from the 
user parameter file and the current date. The header record Is written as 
the first record on the new classification file. 

CLASFYG then evaluates the test field using the following steps: 

A line of pixel data Is read and the greenness value calculated TRJXM for 
every pixel on 'the line Is called for each of the pixels. 

The values for Alpha, Beta, t^ and the CHI-squar» are calculated by TRJXM 
for each pixel . 

The line of pixels coefficients for Alpha, Beta and t^ and the CHI-square 
value Is written to the coefficients file. 

Flowcha rt 
Not applicable. 


Listing 

See Appendix D for program. 
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3.5.2 SUBROUTINE CROPG 


Purpose 

CROPG reads and analyzes cards describing the Initial guesses for A, alpha, 

beta and t^ for the training field, the name of the crop and optionally the 
sun angle. 

Linkages 

CROPG Is called by CLASFYG. 

Interface 

Calling sequence: 

CALL CROPG (lERROR). 

Calling sequence parameters: 

Input/ 

Argument Output De icriptlon 

lERROR 0 0 If no error occurs; 

occurs . 


blocks are contained In Appendix. 
Input/ 

Output Description 

0 Name of crop of Interest. 

0 Initial guesses for A, 

Alpha, Beta and t^. 

None.. 


1 If error 

Function value: 

Not applicable. 

Labeled COMMON parameters: 

Full description of labeled COMMON 

Element 

Label Variable Posi tion 

/USER/ ICR0P(4) 1 

/USER/ ZX(4) 2 

Blank COMMON parameters: 
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Inputs 

Unit Type 

21 Seq. data 

Ou tputs 
None. 

Storage requirement 
Not applicable. 

Description 

Each record In the user defined parameter file Is read and checked to have a 
valid control word In characters 1 through 4 defining the type of Information 
In characters 11 through 80. If the. control word Is Invalid, an error 
message Is written and the error Indicator Is set nonzero; otherwise, the 
Information Is saved In the appropriate common block element. 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 


Description 

User Information file (See Section 3.4.1) 


3-15 


3.5.3 FUNCTION MOVSYM 


Purpos< 

MOVSYM nov«s characttrs from array CARD to a series of left-justi fled words 
in array ITEMP. 

Linkages 

MOVSYM is called by CR0P6. 

Interface 

Calling sequence: 

KOUNT • fiOVSYM (CARD. ITEMP). 

Calling sequence parameters: 

Input/ 

Argument Output Description 

CARO I 80 characters of user supplied data. 

ITEMP 0 Array of left-justi fled words. 

Function Value: 

Name Description 

KOUNT Count of symbols. 

Labeled COMMON parameters: 

None. 

Blank COFflON parameters: 

None. 

Inputs 

None. 


Outputs 

None. 
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Storagt ’‘gqutrcmtnt 
Not applicable. 


Description 

MOVSYM separates the characters In array CARD starting at character 11 and 
creates a series of left-justified words In array ITEMP. 

FI owchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.4 FUNCTION NUMBFL 
Purpose 

NUMBFL changes floating point numbers from characters In array CARD to 
computational numbers In arrary XNMVEC. 

Linkages 

NUM< . Is called by CROPG. 

Interface 

Calling sequence: 

KOUNT • NUWFL (CARO, XNMVEC). 

Calling sequence parameters: 

Input/ 

Argument Output 

CARO I 

XNMVEC 0 

Function value: 

Name Description 

KOUNT Count of floating point numbers. 

Labeled COttlON parameters: 

None. 

Blank COFfON parameters: 

None. 

Inputs 
None. 

Outputs 
None. 
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Description 

80 characters of usfr supplied data. 

Array of floating point computational numbers. 


Storaqt rtgulrtwunt 
Not appllcabit. 


Dtscriptlon 

NUWFL stparatas tha characters In array CARO starting at character 11 and 
creates a series of computational floating point numbers. 

Flowchart 
Not applicable. 

Listing 


See Appendix 0 for program. 


3.5.S ROEXEC 
Purpost 

ROEXEC reads the segment number, number of acquisitions, acquisition dates 
and Landsat numbers. 

Linkages 

ROEXEC Is called by CLASFYG. 

Interface 

Calling sequence: 

CALL ROEXEC. 

Calling sequence parameters: 

None. 

Function value: 

Not applicable. 

Labeled COMMON parameters: 

Full description of labeled COMMON blocks are contained In Appendix A. 


labgi 

Variable 

Element 

Position 

Input/ 

Output 

Oescriptlon 

/USER/ 

SUNANG(8) 

3 

0 

Sun angle for each 

/MI SC/ 

NOACQ 

1 

0 

acquisition. 

Number of acquisitions. 

/MI sc/ 

ACQ0T(8) 

2 

0 

Oate for each acqu1s1t1( 

/MISC/ 

KSEGM 

3 

0 

Segment number. 

/MI SC/ 

LAN0ST(8) 

c 

0 

Landsat number for each 

Blank COMMON parameters: 
None. 



acquisition. 
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Unit Type 


Description 

7 Seq. data Data from RT&£ data base (See Section 3.4.3). 

Outputs 
None . 

Storage requirement 
Not applicable. 

Description 

RDEXEC reads the segment number, number of acquisitions, acquisition dates 
and Landsat numbers. 

Flowcha rt 
Not applicable. 

Listing 

See Appendix 0 for listing. 




3.5.6 SUBROUTINE ROFLES 
Purpose 

ROFLES reads 1) the header record and one set of vertices or 2) one line of 
radiance values. 


Linkages 

ROaES is called by CLASFYT. 


Interface 
Calling sequence: 

CALL ROFLES (I TYPE. lOATA. X. Y, I ERROR). 


Calling sequence numbers: 


Argument 

Input/ 

Output 

Oescription 

. I TYPE 

I 

1 to read header from each acquisition 
file and field definition; 2 to read 1 
line of radiance values from each 
acquisition file. 

I0ATA(856,8) 

0 

Line of radiance values from each 
acquisition file. 

X{4) 

10 

Array of sample vertices defining field. 

Y(4) 

10 

Array of line vertices defining field. 

I ERROR 

10 

0 - no error; 

1 - error. 

Function value: 



Not applicable. 




) 




Labeled COMMON parameters: 

Full description of labeled COMMON blocks are contained in Appendix A. 
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Label Variable 

/LINES/ NOACQS 

/LINES/ NXTLNE 

Blank COMMON parameters: 
None. 

Inputs 


Unit 

Type 


11 

Seq . 

data 

>. 12 

Seq. 

data 

13 

Seq. 

data 

14 

Seq . 

data 

15 

Seq. 

da ta 

16 

Seq. 

data 

17 

Seq. 

data 

18 

Seq. 

data 

Outputs 




None. 

Storage requirement 
Not applicable. 


Element Input/ 

Position Output 

4 I 

1 10 


Description 

Number of acquisitions. 
Next line number. 


Description 

Acquisition data 

Acquisition data 

Acquisition data 

Acquisition data 

Acquisition data 

Acquisition data 
Optional . 

Acquisition data 
Optional . 

Acquisition data 
Optional . 


for date 1 (See 
for date 2 (See 
for date 3 (See 
for date 4 (See 
for date 5 (See 
for date 6 (See 

for date 7 (See 

for date 8 (See 


Section 4.3.2) . 
Sectljn 4.3.2). 
Section 4.3.2) . 
Section 4.3.2) . 
Section 4.3.2) . 
Section 4.3.2) . 

Section 4.3.2) . 

Section 4.3.2). 


Description 

If ITYPE«1, RDFLES reads the reader record from each of the universal 
formatted acquisition data tapes and the vertices for one field. If ITYPE-2, 
RDFLES reads one line of radiance values for each acquisition. 
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Flowchart 
Not applicable 


Listing 

See Appendix 0 


3.5.7 SUBROUTINE RDVERT 
Purpose 

RDVERT reads one set of vertices. 


Linkages 

RDVERT Is ca'lled' by CLflSFYG."' 

Interface 

Calling sequence: 

CALL RDVERT (XFLD, YFLD, lERROR). 


Calling sequence parameters: 

Input/ 

Argument Output 


XFLD(4) 0 

YFLD(4) 0 

lERROR 0 


Description 

Array of sample vertices defining field. 
Array of line vertices defining field. 

0 - no errors; 

1 - error. 


Function value: 

Not applicable. 

Labeled COMMON parameters: 
None. 


Blank C0^W)N parameters: 
None. 
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Inputs 

Unit Type Description 

21 Seq. data User specifUd data file (See Section 4.3.1). 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

RDVERT reads the vertices for one field. 


Flowchart 
Not applicable. 

Listing 

See Appendix D. 


3-26 


3.5.8 SUBROUTINE RECT 
Purpose 

RECT determines if a pixel lies within the specified vertices of a specific 
area. 


Linkages 

RECT is called by CLASFYT. 


Interface 
Calling sequence: 

CALL RECT (LINE, PIXEL, YY, XX, ICOR). 


Calling sequence parameters: 

Input/ 

Argument Output 


LINE I 

PIXEL I 

YY(4) I 

XX(4) I 


Description 

Current line number of pixel. 

Current sample number of pixel. 

Array of sample vertices defining field. 
Array of line vertices defining field. 


ICOR 0 0 if pixel in field; 

1 it pixel not in field. 

Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 
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Inputs 


None. 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

The pixel Is checked to be a boundary pixel and considered In the training 
field If It Is a boundary pixel. If the pixel Is within the minimum and 
maximum limits of the specified area, the angles created between Mnes 
drawn between the pixel and boundary pixels are calculated and summed. 

If the sum of the angles Is 360; the pixel Is considered to be In the 
specified field. 

Flowchart 
Not applicable. 


Listing 

See Appendix D for program. 


3.5.9 SUBROUTINE SOEV 
Purpose 

SOEV calculates the mean and standard deviation of a vector. 


Linkages 

SOEV Is called by CLASFYT. 

Interface 

Calling sequence: 

CALL SOEV (X. N. XB. S) . 


Calling sequence parameters: 


Argument 

X(N) 


N 

XB 

S 


Input/ 

Output 

I 

I 

9 

9 


Oescriptlon 
Array of values. 

Number of values In array X. 

Mean of values In array X. 

Standard deviation of values In array 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 


Inputs 

None. 
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Outputs 

None. 


Storage requirement 
Not applicable. 

Description 

SDEV calculates the mean (XB) and standard deviation (S) of vector X with 
N points. 

Flowchar t 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.10 SUaPOUTINE TRJXM 
Purpose 


TRJXM makes a least-tjquare fit to a nor-1 Inear function to calculate 
coefficients for the curve. 


Linkages 

TRJXM Is called by CLASFYG. 


Interface 
Calling sequence: 

CALL TRJXM (X, Y, SIGMAY. NFILES, NTERMS, MODF, A DELTAA, SIGMAA, 
FLAMOA, YFIT, CHISQR). 


Calling sequence parameters: 

Input/ 

Argument Output 


X I 

Y I 

SIGMAY I 

NFILES I 

NTERMS I 

MODE I 

ESTIM I.O 

DELTAA I 

SIGMAA 0 

FLAMOA 0 

YFIT 0 

CHISI,R 0 


Descriptio n • 

Array of Julian day/100 for acquisitions. 
Array of dependent variables. 

Array of errors on Y. 

Number of acquisitions. 

Number of degrees of freedom. 

Not used, was switch for calculation weight. 
Array of coefficients A, Alpha, Beta and t^. 
Array of Increments for A. 

Standard deviations. 

Proportion of gradient search. 

Array of fitted values. 

Residual CHI square value. 
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Input/ 


Argument 

Output 

Description 

IND 

N/A 

Not used. 

WEIGHT 

I 

Array of weighting constants. 

IPRINT 

I 

Print flag to use In debug. 


Function valut: 

Not applicable. 

Labeled COftVN parameters: 

None. 

Blank COMMON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

TRJXM calculates alpha and beta matrices, using the least-square method, 
for defining all parameters of the curve for the selected pixels In the 
training field. 

The CHI-squarc value of the starting point Is evaluated and the modified 
curvature matrix Is Inverted to find the new parameters. 

The convergence depends on the CHI-square value and on the size of step. 
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If the CHI-square value Increases, FLAMOA Is Increased and the process 
Is re-lnitlated. If the CHI>square value descreases, the process Is 
re- Initiated until convergence Is obtained. The parameters and uncer- 
tainties are evaluated. 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.11 SUBROUTINE F2ERVG 
Purpose 

F2ERV6 computes the derivatives of the function of each of the X values. 
Linkages 

F2ERVG Is called by TRJXN. 

Interface 
Calling sequence: 

CALL F2ERVG (X. I, ESTIM. DELTAA, NTERMS, DERIV) . 


Calling sequence parameters: 

Input/ 

Argument Output 


Description 


X 

I 

ESTIM 

OELTAA 

NTERMS 

DERIV 


I Array of time points for acquisition. 

I Current Index In calling program. 

I|9 Array of coefficients A, Alpha, Bet? 

and t^. 

I Array of Increments for A. 

I Number of degrees of freedom. 

0 Array of derivatives with respect to 

parameters ESTIM. 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 
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blank CCtflON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 

Storage requiremen t 
Not applicable. 

D escription 

The radiance value and am set positive. Y Is calculated as the 
radiance value/t^. If Y Is less than or equal to 1, the 4 derivatives 
are set to very small numbers; otherwise, the derivatives are set to 
calculated values. 

Flowchart 
Not applicable. 

Listing 

See Appendix 0 for program. 
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3.5.12 FUNCTION FCHISQ 
Purpose 

FCHISQ calculates the CHI-square value per degree of freedom. 
Linkages 

FCHISQ Is called by TRJFIT and CHNFIT. 

Interface 
CalUng sequence: 

GOODFT-FCHISn (Y. SIGMAY, NPTS, NFRCF, MODE, YFIT, WEIGHT). 
Calling sequence parameters: 


Input/ 


Argument 

Output 

Description 

Y(NFILES) 

I 

Array of dependent variables. 

SIGMAY(NFILES) 

I 

Array of errors on Y. 

NPTS 

I 

Number of acquisitions. 

NFREE 

I 

Number of degrees of freedom. 

MODE 

I 

Not used. 

YFIT(NFILES) 

I 

Array of fitted values. 

WEIGHT(NFILES) 

I 

Array of weighting factors. 

Function value: 



Name Description 



GOODFT CHI-square value per degree of freedom. 

Labeled COMi'lON parameters: 

None. 


Blank COMMON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

FCHISQ calculates CHI-square value - Ki YFIT-Y)/SIGMA)^ 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.13 FUNCTION F2Ni;TG 


Purpose 

F2NCTG calculates tne value of DERIV from the radiance value, alpha and 
beta . 

Linkages 

F2ERVG Is called by TRJXM. 

Interface 

Calling sequence: 

CALL F2NCTG (X, I, ESTIM, IND) . 

Calling sequence parameters: 

Input/ 

Argument Output 
X(8) I 

I 1 

ESTIM(4) I 

IND N/A 

DELTAA(8) I 

NTERMS I 

0ERIV(4) 0 

Function value: 

Name Description 

YFIT(8) Array of fitted values. 

Labeled COMMON parameters: 

None. 


Description 

Time array of Julian date/100 for each 
acquisition. 

Index In calling program. 

Array of coefficients A, Alpha, Beta and t^. 
Not used. 

Array Incremental values. 

Number of degrees of freedom. 

Derivatives . 
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Blank COMMON parameters: 

None. 

Inputs 
None . 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

F2NCTG calculates the value of the points on the curve from the values of 
alpha and beta and the current radiance value. 

Flowcha rt 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.14 SUBROUTINE MTINVl 


Purpose 

^tTINV^ Inverts a matrix of order NORDER. 

LI nkaqes 

MTINVl is called by TRJXM. 

Interface 

Calling sequence: 

CALL MTINVl (ARRAY, NORDER, DET) . 

Description 
Matrix. 

Order of matrix. 
Determinate of matrix. 

Function value: 

None. 

Labeled COMMON parameters: 

None. 

Blank COMMON parameters: 

None. 

Inputs 

None. 


Calling sequence parameters: 


Input/ 

Argument Output 

ARRAY(NORDER) 10 

NORDER I 

DET 0 


Outputs 

None. 
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storage requirement 
Not applicable. 


Description 

MTINVl uses the Gauss-Jordon reduction techniques. 

Flowchart 
Not applicable. 


Listing 

See Appendix D for program. 


3.5.15 FUNCTION NUMBR 
Purpose 

NUMBR changes Integers from character array CARO to computational numbers In 
array NUMVEC. 

Linkages 

NUMBR Is called by CROPTM. 

Interface 
Calling sequence: 

KOUNT - NUMBR (CARD, NUMVEC). 


Calling sequence parameters: 


Argument 

Input/ 

Output 

Description 

CARO 

I 

80 characters of user supplied data In character 
format. 

NUMVEC 

0 

Array of Integer computational numbers. 

Function value: 



Name 

Description 


KOUNT 

Number of computational Integers 

Labeled COMMON parameters: 



None. 

Blank COMMON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 
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storage requirement 
Mot applicable. 


Description 

NUMER separates the characters In array card starting at the 11th character 
and creates a series of computational Integers. 

Flowchart 
Not applicable. 

Listing 

See Appendix 0 for program. 
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4.0 OPERATION 


4.1 OPERATING DESCRIPTION 

CLASFYG Is operational on the IBM 3031 computer at '.ARS, West Lafayette, 
Indiana. 

The CLASFYG program Is one of the programs of the BADHWAR SYSTEM. 

CLASFYG requires the use of a 0 disk which Is assigned as a temporary disk 
and an E disk which Is used to load LARS routines onto. The user therefore, 
must not assign a disk to his machine using either MODE E or MODE 0. These 
disks will be assigned as needed. 

Prior to executing the CLASFYG program the user must establish on his A disk 
a USER INFORMATION file as described In Section 3.4.1. 


4.2 COMMANDS DESCRIPTION 


To execute CLASFY6 the user will enter the following series of comnands 
which Invoke the JOB CONTROL SOrTWARE. These commands are divided Into two 
classes namely (1) FUNCTION commands and (2) PROGRAM commands. The FUNCTION 
commands* which perform all the functions except executing the program, are 
reusable, I.e., once they are Invoked they remain In effect until reentered. 
The PROGRAM comnands, which simply execute the program, must be reentered 
each time the program Is to be executed. 

The following list gives the commands required to execute the CLASFYG 
program. They an all FUNCTION commands except the PROGRAM command 
CLASFYG. These commands are to be given In the listed order. 

START 

CLASFYG 

END 

The following sections describe each of the commands In detail. Input 
fields are separated by blanks. If more than one word Is required to 
describe an Input field, the description Is enclosed In pointed brackets <>. 
If an input Is optional the field Is enclosed In square brackets []. Do 
not Include these explanatory characters <> [] when actually submitting Input 
to the computer. To enter a command the user types one Input per defined 
Input field and separates each field with a blank. 


4.2.1 START 


The START command spools the user's console file. The use of this comnand 
along with the END command will provide a listing of all Information 
appearing on the user console file. (If running an Interactive job this 
Is the terminal.) If running a batch job this Is a system defined device. 
The START command Is Invoked by the user typing the following; 

START 
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4.2.2 CLASFYG 

Tht CLASFYG connand Is a PROGRAM command Is used to Invoke the CLASFYG 
program. All landsat Input data Is from the LARS RT&E Data Base and a 
series of programs Is required to perform the necessary Interface. The 
following Illustrates this software flow. 

..RTEER (LARS ROUTINE) 

CLASFYG LCGINF i.SEGALL (LARS ROUTINE) 

!.TSRTHS....IVALUE 

For a detailed descriptloii of the above JOB CONTROL SOFTWARE see appendix 
C. 

The CLASFYG connand must not be given unless the user has established a 
user's file on his A disk under the file description <SEGHENT NUMBCR> USERG 
as described In Section 3.4.1. Furthermore.the SEGMENT NUMBER must contain 
all 4 digits of the segment. For example: segn^^^lt 0882 dictates that the 

USER INFORMATION file be named: 

0882 USERG A 

The CLASFYG connand Is Invoked by the user typing the following. 

CLASFYG <SEGMENT> <ACQ.#1> <ACQ.#2> <ACQ.#3> <ACQ.#4> <ACQ.#5>. . .<ACQ.#8> 

As described In Section 3.4.2, the CLASFYG program requires the user to 
specify a minimum of 5 segment acquisitions as Input and a maximum of 8 
may be specified. 

The classification file output from an execution of CLASFYG Is written to 
a file named <SEGMENT NUMBER> COEF A. The program Inputs acquired from 
the LARS Data Base and the USER INFORMATION file Is spooled to the HOUSTON 
printer. 
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4.2.3 END 


This conmand closts the user's console file and causes a spooled copy to be 
sent to the HOUSTON printer. This conmand has no effect If the START 
conmand was not previously Issued. The END conmand Is Invoked by the user 
typing the following: 

END 
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4.3 OPERATING EXAff>L£ 

For our example we assume that the user has established on his A disk the 
required USER INFORMATION file. In addition since our example deal; with 
an execution using segment 882 the USER INFORMATION file must be established 
under the filename 0882 and flletype USERG. 

COMMAND EXPLANATION OR ACTION TAKEN 


START Spools the console file. 

CLASFYG 0882 78150 78159 78186 78313 78221 78231 78267 

Executes the CLASFYG 
using 7 acquisitions of 
segment 882. 

END Closes the user's console 

file and spools the file 
tothe HOUSTON printer. 
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APPENDIX A 


COMMON BLOCK /LINES/ 

NXTLNE Not currently used. 

LINENO Current line number. 

IDATA(856,8) Acquisition data. 

NOACQS Number of acquisitions. 
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COMMON BLOCK /MISC/ 


NOACQ 

ACQDT(8) 

KSE6M 

LANDST(8) 


Number of acquisitions* 

Acquisition dates stored left Justified In double precision 
word. 

Segment number stored as 4 character word. 

Landsat numbers. 
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COffICN BLOCK /USER/ 


ICR0P(4) 

ZX(4,10) 

SANGGG(8) 


Name of crop of interest. 

Initial guesses for A, Alpha, Beta and t^ for greenness. 
Sun angles for acquisitions. 
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APPENDIX B 


CLASFYG Variable Definitions 


Variable 

i='BS 

FiCPDT 

flCODTE 

t=iLPHEF 

PNG 

BETPBF 

BtEG 

CHIBUF 

CHISO 

CHNVRL 

CLNPST 

CPTINE 

DBUF 

HDP 

IftLPHR 

IBETP 

ICHISO 

ICDP 

ICPTIM 

ICPDP 

ICUT 

IDPTR 

IPPTE 

IDPY 

IDIFFC 

IDIFFV 

TDUPDT 

IHDLPl 

IH0LD4 

II 

ILIM 

I MON 

INCPE 

INDEX 

IPIXEL 

ITD 

ITYPE 

IVCPTN 

lYERP 


Definition 

PBSDLUTE VPLUE FUNCTION 

RCOUISITION DRIES IN RLPHRNUNEPIC CHRPRCTEP 
RCOUISITION DRIES IN 1 CHRPRCTEP WOPDS 
BUFFEP FOP RLPHR VRLUES 

TEMP VRLUE. USED IN SUN RNGLE COPPECTION 

BUFFEP FOP BETR VRLUES RFTEP TEST PIXEL CRLL TO TPJXM 

SEGMENT NUMBEP IN CHRPRCTEPS 

BUFFEP FOP CH SO. VRLUES 

FP CHI SO. VRLUE CRLC. BY TPJXM FOP CUPPENT TEST PIXEL 

HOLD RPER FOP CHRNNEL PRDIRNCE VRLUES 

CONSTRNT FOP CONVEPTING LRNDSRT 3 DRTR 

SUBROUTINE TO PETUPN TODRY^S DRTE RND TIME 

NOT USED 

BUFFEP FOP HERDEP 

2ND RPPPX VRLUE* RLPHR. CRLC BY TPJMX FOP TEST PIXEL 

3PD RPPPX VRLUE* BETR. CRLC BY TPJMX FOP CUPPENT TEST PIXEL 

INTGP CHI SO. VRLUE CRLC. FOP CUPPENT TEST PIXEL BY TPJXM 

RETURN PRPRMETEP FROM PECT* 0-PIXEL IN FLD* l-PIXEL NOT IN FLD 

INITIRL CP TIME IN HUNDREDTHS OF SECONDS 

CROP IN RLPHRNUMEPIC CHRPRCTEPS 

NOT USED 

BUFFEP RPER FOP LINES OF PRDIRNCE VRLUES 
RCOUISITION DRIES IN INTEGER FOPMRT 
DRY IN CHRPRCTEP FOPMRT 

ELRPSED CP TIME IN HUNDREDTHS OF SECONDS 
ELRPSED VIPTURL CP TIME IN HUNDREDTHS OF SECONDS 
INDEX TO DUPLICRTE DRTE. ZERO IF NO DUPLICRTES 
RPPRY OF 1 CHRPRCTEP WOPDS 
RPPRY OF 4 CHRPRCTEP WORDS 
INDEX 

SET* BUT NOT USED 

MONTH IN CHRPRCTEP FOPMRT 

TEMP VRLUE USED IN UNPRCKING PIXELS 

INDEX TO LRNDSRT COPPECTION FRCTOP 

CUPPENT PIXEL 

4TH RPPPX VRLUE* TO* CRLC FOP CUPPENT TEST PIXEL BY TPJXM 
INPUT TO SUBP. PDFLE TO PERD HERDEP a> OP PRDIRNCE VRLUES»2) 
INITIRL VIPTURL CP TIME 
2 CHRPRCTEP YERP 
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VARIABLE DEFINITION 


JM 

JtEGM 

► BB 
••'BLFiMK 
>■ EOL'PL 

► DPr<EP 
K"iEi5M 

► IJPRCO 
KUPCHM 
KUPPXL 
LR 

LRMPST 

LIMIT 

LI NEMO 

LMRX 

LMIh 

LDC 

MRXPXL 

Mir<PXL 

MODE 

NCPTIM 

MMGDE 

NORCO 

NDRCOS 

NORCOS 

NOPXLS 

NOTPPy 

NTEPM 

NT I ME 

NVCPTM 

NXTLNE 

PELTR 

OELTR 

OFTX 

PDFLES 

PECT 

SANG 

SRNGGG 

SDEV 

SETIJPG 

SGMRTP 

SIGMR 

SIGMRG 

SIN 

ST I ME 


original 
OF POOR 

INDEX 

SEGMENT MJMBEP RS 1 CHRPRCTEP l.iOPDS 
INDEX 

CONSTANT OF BLANK 

CONSTANT OF EOURL MRPK 

INDEX TO VERTICES FOP TEST FIELD 

SEGMENT NUMBER AS 4 CHRPRCTEP UDPDS 

CUPPENT ACOUISITION 

CURRENT CHANNEL 

CUPPENT PIXEL 

INDEX 

RPPAY OF LANDS AT NUMBERS FOP ACOUISITIONS 

MAX CALLS TO TPJXM FOP EACH PIXEL 

CUPPENT LINE NUMBEP 

MAXIMUM LINE NUMBEP FOP BOTH FIELDS 

MINIMUM LINE NUMBEP FOP BOTH FIELDS 

TEMP VALUE USED IN UNPACKING RADIANCE VALUES 

MAXIMUM PIXEL 

MINIMUM PIXEL 

NOT USED. ORIGINALLY USED TO SELECT WEIGHTING METHOD 

ENDING CP TIME IN HUNDREDTHS OF SECONDS 

CONSTANT OF 0. USED IN CALL TO TPJXM 

NUMBEP OF ACQUISITIONS 

NUMBEP OF ACQUISITIONS 

NUMBEP OF ACQUISITIONS 

NUMBEP OF PIXELS 

CUPPENT NUMBEP OF PIXELS IN TPAINING FIELD 
NUMBEP OF DEGPEES OF FREEDOM <4> 

CUPPENT DATE. NOT USED 

ENDING VIRTUAL CP TIME IN HUNDREDTHS OF SECONDS 
NOT USED 

8TH PARAMETER TPJXM. CRN BE USED IN NON-RNRLYTIC DIFFERENTIATION 
NOT USED 

IITH PARAMETER IN CALLING SEQUENCE TO TPJXM 
SUBROUTINE TO READ FILE DATA 

SUBROUTINE TO DETERMINE IF PIXEL IN USER DEFINED FIELD 
SUN ANGLE ARRAY 
USER INPUT SUN ANGLE 

SUBROUTINE TO CALCULATE MEAN AND ST. DEV. 

SBUPOUTINE TO READ USER INPUT DATA 

ST. DEV. GPEEN TRAIN FLD. INDEXED BY ACQ. NO. USED FOP TEST FLD 
ERROR ON PARAMETERS FOR RETURN FROM TPJXM 
ST. DEV. FOP TPAINING FIELD 
SINE FUNCTION 

STARTING TIME OF DAY IN CHARACTERS FOP HEADER 


PAGE IS 
QUAUTY 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


k 


VARIABLE DEFINITION 

TIMECP ELAPSED CP TIME Ih MINUTES 

TIME VP ELAPSED VIRTUAL CP TIME IN MINUTES 

TODAY TODAY'S DATE 

TPAINA 1ST APPPX VALUE CALC IN TPJXM FQP TRAIN FLD. USED IN TE'T TPJXM 

TPF.LPH 2ND APPPX VALUE CALC. IN TPJXM FDR TRAIN FLD« USED IN TEST TRJl'M 
TRBETA SRD APPPX VALUE CALC IN TPJXM FDR TRAIN FLD* USED IN TEST TPJXM 

TPJXM SUBROUTINE TO CALCULATE FIT ON TRAINir<i3 FLU AND THEN TEST PIXEL 

TPTO 4TH APPPX VALUE CALC IN TPJXM FQP TRAIN FLD* USED IN TEST TPJXM 

TZEPO 4TH APPkX VALUE. TO. CALC FOP CURRENT TEST PIXEL BY TPJXM 

TOBUF BUFFER FDR TO ARRAY AFTER TEST PIXEL CALLS TO TF JXM 

VEPTPX USER SUPPLIED LINE VERTICES 

MEIGHT 1./ SH5MAY»»2f WEIGHT FACTOR USED BY STATISTICAL ROUTINES 
>:BETA FP 2ND APPPX VALUE CALC. BY TRJMX FOP CURRENT TEST PIXEL 
XCHIN INITIALLY 200. THEN CURRENT VALUE OF CHI SQ AFTER CALLS TO TPJXM 

XDAY FL. PT. VALUES FOP DAYS OF YEAR / 100. 

XDIF FL. PT. TEMP HOLD AREA 

XGREEN GREENNESS VALUES FOR TRAIN PIXELS 

XI SOP CHI SO. AFTER BOTHS CALLS TO TPJXM 

XLAMDA INCPE PARAMETER. SET l.E-03 1ST CALL TO TPJXM. MODIFIED BY TPJJM 

XLINE FP USER SUPPLIED LINE NUMBER FDR CALL TO PECT 

XMEAN MEAN OF GREENNESS VALUES FDR TRAINING FIELD 

XMNTPN ARRAY OF MEANS FDR TRAINING FIELD ACOUISITIONS 

XPIXEL FP USER SUPPPLIED PIXEL NUMBER FOP CALL TO RECT 

XPIXGP ARRAY OF GREENNESS VALUES FDR PIXELS 

XTH TEMP 

XTOL TOLERANCE FACTOR USED FDR TPJXM FIT CHECK 

VEPTLN USER SUPPLIED PIXEL VERTICES 

VGPEEN GREENNESS VALUES FDR TRAINING FIELD 


L 
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